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We study an exactly solvable model describing a stripe consisting of a Toda array of N anharmonic 
elastic chains sandwiched between two conducting chains. It is shown that the presence of a charge 
on one chain generates a gapless excitation branch (Luttinger liquid) on the other and leads to 
increase in the phonon frequencies. 
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During recent several years there have been rather in- 
tense interest in systems which combine one-dimensional 
and two- or three-dimensional features. As examples of 
such systems one can mention domain wall (or stripe) 
formations in doped Mott insulators jjj and quasi-one- 
dimensional systems with electron-phonon interactions. 
The capital fact about phonon systems is that phonons 
are never one-dimensional. Thus in both cases electrons 
residing on one-dimensional chains are immersed in a 
two- or three-dimensional active environment. This fea- 
ture makes such problems very challenging from the the- 
oretical point of view because on one hand the low di- 
mensionality leads to the enhancement of interactions 
thus requiring non-perturbative approach, and on the 
other hand the multi-chain structure fits rather awk- 
wardly with available non-perturbative techniques. Usu- 
ally theorists resort to a hybrid treatment incorporating 
features of exact solution for individual chains and RPA 
(Random Phase Approximation) for interchain interac- 
tions 

From the qualitative point of view there are several 
interesting problems concerning the systems in question. 
One is how conducting regions affect each other and an- 
other is how changes in the conducting region affect the 
active medium in which it exists. 

Here we discuss these problems using as an example 
an exactly solvable model of a static stripe suggested 
by Fateev j|. A remarkable property of this model is 
that it exists in two different representations describ- 
ing situations with different physics. One representa- 
tion is fermionic and the other one is bosonic and they 
are related by a duality transformation. In both its 
incarnations the model describes two species of (1+1)- 
dimensional charged particles located on different sides of 
an insulating stripe and interacting via its elastic modes 
(see Fig. 0). The latter ones are anharmonic optical 
phonons described by the so-called Toda array. The La- 
grangian density has the following form: 



N) = Croda + C-edge 
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In its fermionic incarnation the edge Lagrangian is given 
by 



42e= E [^.7^. + Y^'^.) a 



m^2e^ N - {m/(3f (e" Wl +e 2Ww ) 



(3) 



where g = P 2 /(Att + 1 ). It can be interpreted as a 
model of two charge density waves coupled through the 
Toda array. 
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FIG. 1. The system under consideration: two conducting 
regions coupled by the Toda array. 

The bosonic representation is more intriguing because 
of its connotations with the hypothetic stripe formations 
in doped Mott insulators. The edge Lagrangian density 
is 
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e 2^1 + (/3/2) 2 A*A s 



^-(|Ax| 2 + 2)e" Wl - ^-(|A 2 | 2 + 2)e^~ (4) 

where A s are complex bosonic fields which can be inter- 
preted as superconducting order parameters residing on 
the edges of the stripe. 

The duality transformation relating fermionic to the 
bosonic representation is 



1 



L^((3,N) =L^{A-k/P,N -2) 



(5) 



which also corresponds to g — » 1 — g. Therefore the 
strong coupling limit of one model is related to the weak 
coupling limit of another one. 

The model has U(l)xU(l) symmetry and the corre- 
sponding charges are given by 



dx 



(A* s doA s - A s d A* s 
[l + (/?/2)2|A s |2] 



(6) 



One can modify the Hamiltonian introducing chemical 
potentials coupled to the charges H — > H — h + Q + — 
h-Q^ which would correspond to populating the edges 
of the stripe. In this publication we shall be particularly 
interested in the situation when one chemical potential 
is much greater than another. 

In H Fateev gave a proof of integrability of the above 
models and suggested exact two-body S'-matrices. At 
zero chemical potentials all excitations have gaps. The 
particles living on the edges may be considered as fun- 
damental ones, excitations of the Toda array being their 
bound states. There are N — 1 bound states with the 



Mj = 2Msm(%j/2X),(j 



l,2,...N-l),X = N-g (7) 



where M = M(m, (3) is the mass of the fundamental 
particle. All bound states are neutral. The case N = 1 
where there are no bound states has been thoroughly 
studied in || . 

Each fundamental particle is labeled by two indices 
corresponding to the two U(l) groups (it is either soliton 
or antisoliton and it belongs to one chain or another). 
The two-body scattering matrix of the fundamental par- 
ticles is given by a tensor product of two sine-Gordon 
S'-matrices multiplied by a CDD-factor responsible for 
calculation of double poles. Since the Bethe ansatz equa- 
tions for the sine-Gordon model is well known, it is quite 
straightforward to write down such equations for the 
model in question: 
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(10) 


i=i 


6=1 





where e n (x) = sinh X(x + inir/2)/ sinhA(x — mm/2), the 
numbers ra\^. = n/2 — Qi^ and 



N-l 



q (a\ 2i5 A (e) TT sinhfl - isin(7ra/A) 

S Q (d)=e M ; I I -.a,--/ 7TT Wlth 
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6 X (0)= / dc^ 
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sm(w6) sinh[7rw(A- 1 - l)/2] 
oj cosh(7rw/2) sinh(7ra;/2A) 



We claim that that creating a charge on one shore of 
the stripe one generates a band of soft collective modes 
on the other shore. This statement can be understood 
without detailed calculations which we leave for the ex- 
tended publication, just using a qualitative analysis of 
the Bethe ansatz equations (|| p|jic|) and will be illus- 
trated on figure ||. As a matter of fact, similar effect may 
occur in all integrable systems where the S-matrix is a 
tensor product (here Principal Chiral Fields may serve 
as a good example - see, for example, [||). In all these 
cases the chemical potential coupled to one of the con- 
served charges generates not only its own charge, but also 
induces a finite susceptibility with respect to the other 
chemical potential. Thus the second subsystem becomes 
gapless. We call this effect the induced Luttinger liquid. 

Notice that Eqs. (p|JlO|) for u a , v a look very similar 
to BA equations for the spin- 1/2 XXZ spin chain (see, 
for example, ||) with the anisotropy related to the pa- 
rameter g. In the presence of a large chemical potential 
h- <C M applied on the left shore, the density of particles 
in the Fermi sea n/L is finite and the real rapidities 6 in 
Eqs.(p|[9|]l0|) with energie E are densely distributed over 
some finite interval [-B, B]. The excitation energies e_„ 
of the u-particles are pushed up by the chemical poten- 
tial. The excitation energies e n of the v- particles have 
their minima at \v\ — > oo and therefore at T <c h- (the 
Fermi energy of the ^-particles) they are located very far 
away in the momentum space from the Fermi momenta 
of the ^-particles. As all energies e_„ are gapped,the ra- 
pidities u drop out of consideration and in the equations 
for v one can neglect a feedback between v's and #'s and 
replace 6's by their ground state values. Furthermore, at 
large values of \v\ <C B which are important at low ener- 
gies in these equations one can neglect #'s in comparison 
with v's. Then Eqs. (|To|) become exactly like the Bethe 
ansatz equations for the spin-1/2 XXZ spin chain with 
the total number of sites n equal to the number of 6's in 
the ground state: 
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H = JY}S*S X J+1 + S]S] +1 + cos(ng)Sj Sj+i] (11) 
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where we set < g < 1 using the fact that the region g < 
is mapped onto g > by the transformation 1—g — > —g. 
Since all these arguments work only for small energies 
where the spin chain is equivalent to the Luttinger liquid, 
one can say that the ^-subsystem describes a spinless 
Luttinger liquid with the Luttinger parameter 



tf = [2(l-s)]- 



(12) 



We emphasise that this value of K is independent of num- 
ber of chains or the number of particles in the Fermi sea. 



2 



FIG. 2. Energies e„(0) before (left) and after (right) application of a chemical potential h- on the left wire. Only those 
modes are shown whose energies do not vanish at T — 0; the chemical potential creates large gaps ~ (n — l)/i_ for (n I 2) 
and smaller gaps ~ A p h (p0[) for the phonon modes Ej (j = 1, ... N -1) and e_i. The mode ei becomes gapless. 



The corresponding velocity V s however, does depend 
on N and the number of particles. One can show that V s 
vanishes at g — ► 0, 1 indicating that the collective band 
described here is generated by the interactions. 

Now we shall briefly outline the technical details lead- 
ing to the above conclusions. Following the standard pro- 
cedure we can derive the Thermodynamic Bethe Ansatz 
(TBA) equations corresponding to Eqs.(pp|JlO|). To make 
the problem simplier we use the fact that the above Bethe 
ansatz equations are invariant under the transformation 



-g=l-g',N = N'-l 

u a X = u a X + m/2, v a X = v a X + in/2 



(13) 



Due to the periodicity of the phase factors in Eqs.([s|,|9||lC|) 
one can make a replacement 

sinh X(u + in7r/2) sinh(Au — mng/2) 
sinh X(u — inn/2) sinh(Au + mng/2) 

These two facts allow us to restrict our consideration by 
the region 1 > \g\ > 1/2. To further simplify our presen- 
tation we shall give explicit expressions only for the case 
g = \/v > with v > 2 being an integer number where 
classification of solutions is simplier. The free energy is 
then given by 

1 N ^ f 

F/L = --yMj / d0cosh01n[l + e~^ (e)/T ] (14) 

27T ~ — ' / 

i=o J 

where M = M and Mj are given by Eq.(^). 
The TBA equations for the bound states are 

T ln[l + e E ^)/ T ] - G jk *Tln[l + e - E *W r ] = Mj coshtf 
+ G j0 *Tln[l + e~ fio(f,)/T ], (j, k = 1, ...N - 1) (15) 

where the star denotes convolution 

/>oo 

f*g(9) = / d6'f(6-8')g(9>) 



In this publication we do not need explicit expressions for 
the kernels Gjk, Goj', the important features of Eqs.(|l5|) 
being that Gjo > and the energies Ej are directly cou- 
pled only to Eq. The latter property is related to the fact 
that bound states are neutral with respect to both U(l) 
groups. 

When the chemical potentials /i_ and h + are applied 
respectively to the left and right edges, the full system of 
TBA is given by 



Tln[l 



] -K*Thx[l 



-E ( 



9)/Tj 



= M cosh 9-(h+ + h-)/2 + G 0j * Tln[l + e"^ {e)/T ] 
- T ^ a n * ln[l + e—" W/ T ] - Ta v ^ * ln[l + e" e -" 

n=l 

-T.s*ln[l + e ei ( 9) / T ] (16) 

equations for e n and e_„ are 

e ±u = h±v/2 — s*T ln[l + c £ ±(-2)/t] ( 17 ) 
e± n = S\ n \,»-ih±v/2 + s * Tln[l + C -^-^ T ][l + e e ±^ T ] 
+ ^-2.MS*Tln[l + e- eWT ] 

-*| n | i i**Tln[l + e- Bo/r ] (18) 
where 

. , sinh[7r(l — n/v)u>/2X} . , r „ , , i 
«nH = \ r ,s(lj) = 2 cosh (iru 2Xv 



sinh[vrcj(l - g)/2X] cosh[vrcj(A + 2g)/2X] 
2 cosh7rw/2 cosh(<77rw/2A) sinh(7ru;/2A) 



The general structure of the TBA can be illustrated by 
the incident diagram drawn on Fig. |j| From these equa- 
tions it is manifest that at h + = and h— > M when Eq 
has a negative part, all e_„ are positive except, possibly 
e_i ( see Fig. f|. 
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FIG. 3. The incident diagram for the TBA equations given 
in the text. The circles correspond to energies and links to 
nonvanishing kernels. The circles on the vertical line corre- 
spond to the bound states. 

From Eq.([l6|) we see that at h- > M the energy Eq 
is negative on some interval [—B,B]. From Eq.(|l5|) it is 
clear that E/s are pushed further up. As follows from 
@ at T = 



E Q w R- 1 * (Mcosb.0 - h-/2) 



(19) 



Substituting this into Eq. ( |15|) we find that the mass gaps 
for bound states are of order of 



A ph ~ M{h-/Mf, v = 9/(N + g) 



(20) 



One can also check that all e n except of e_i still have a 
gap of order of h- and e_i has a much smaller gap ~ A p h- 
Thus at h- 3> M there are two temperature scales in 
the problem: h- itself which play a role of the Fermi 
energy; at temperatures T <C h— one can omit e_„'s from 
Eq. ( |l6|) . The other scale is the much smaller one - it is the 
phonon frequency scale A p h- It also follows from Eqs.jl^) 
that all e n 's with n > 1 are positive and of order of T 
and ei is purely negative in the limit T = 0. Therefore 
the term with ei in the right hand side of Eq.([l6|) gives 
no contribution at T = 0. 

/From now on we shall ignore bound states assum- 
ing that T <C A p h- At such temperatures we can ne- 
glect temperature dependence of B. We are interested 
in the region where all e n are small. This corresponds to 
\6\ — > oo which are far from the Fermi surface of Eq(0). 
Therefore we can replace — Tln[l+e~ E °^^ T ] in the right- 
hand side of Eqs.([l7 18) by its value at zero temperature 



e (-)(6>) = E (6)@[-E (6)]. For 6 > B the driving term 
can be replaced by its asymptotics: 



t-\0)Ms(6) dyE (y)cosh(Xy/g) 



After this substitution, the TBA equations (|lj,|l8j) co- 
incide with the low temperature limit of the TBA equa- 
tions for the anisotropic spin-1/2 Heisenberg magnet (see 
Eq.@ above). 

In order to obtain the free energy at low T we isolate 
contributions which vanish at T = 0. This is achieved 
using the standard TBA machinery; the free energy is 
represented as a sum of two parts corresponding to con- 
tributions from two gapless modes: 



fx = 



h = 



-T I [a^\d)+a { -\e)]\n{l + e-\ Eo ^l T )dd 
-n 2 T 2 /3V c and (21) 



—T J d0s*a-(9) ln[l + e £1 W/ T ] 
~ir 2 T 2 /W s . 



(22) 



This describes independant gapless modes on the left 
edge (/i) as well as the induced mode on the right shore 
(/ia). The function a (8) satisfies 



r(+) 



(0) 



/B yt 
dO'K{0 -0><->(0') = — coshO 
-B 2n 



where cr(+)((9) = at |0| < B and a^(6) = at \6\ > B. 
The latter function is the ground state density of real 0's. 

For very large chemical potential h- ( B — > oo), we 
find the velocities V c = 1 and the velocity of the inducted 
mode V s = g/ (N — g) . Therefore the inducted gapless 
band collapses at g — > as it should. The similar thing 
happens at g — ► 1. 

In conclusion, charging of one edge of the stripe creates 
a band of gapless excitations on the other edge and stiff- 
ens the phonon modes inside. Their energies, however, 
remain much smaller than the Fermi energy which means 
that the system in question has an intermediate energy 
scales where one should expect nontrivial crossovers. 
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